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Abstract 

We compute the two-loop corrections to the thermodynamical pressure of 
an SU(2) Yang-Mills theory being in its electric phase. Our results prove that 
the one-loop evolution of the effective gauge coupling constant is reliable for 
any practical purpose. We thus establish the validity of the picture of almost 
noninteracting thermal quasiparticles in the electric phase. Implications of our 
results for the explanation of the large-angle anomaly in the power spectrum of 
temperature fluctuations in the cosmic microwave background are discussed. 
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1 Introduction 

In one of us has put forward an analytical and nonperturbative approach to 

the thermodynamics of SU(N) Yang-Mills theory. This approach self-consistently 

assumes the 'condensation' of (embedded) SU(2) trivial-holonomy calorons ^ into 

a macroscopically stabilized adjoint Higgs field in the deconfining high-temperature 

phase of the theory ^. This assumption is subject to proof which we establish 

in an analytical way in [H^. The incorporation of nontrivial-holonomy calorons |31 

m ini Hj into the ground-state dynamics can be thermodynamically achieved in an 

exact way in terms of a macroscopic pure-gauge configuration. We thus describe 

the effects of dissociating nontrivial-holonomy calorons (magnetic monopoles which 

attract or repulse one another for small or large holonomy, respectively where 

the former possibility is far more likely.) on the pressure and the energy density 

of the ground state in an average fashion, that is, thermodynamically. By a global 

^2,cioc degeneracy of the ground state and a nonvanishing expectation value of the 

Polyakov loop it can be shown analytically that the electric phase is deconfining. 

Moreover, the inrafred problem of thermal perturbation theory is resolved by a 

nontrivial ground-state structure giving rise to a mass for gauge-field fluctuations 

off the unbroken Cartan subalgebra. 

On tree-level, excitations in the electric phase are either thermal quasiparticles 

or massless 'photons'. The evolution equation for the effective gauge coupling e 
^This phase is referred to as electric phase in jj. 
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in the electric phase is derived from thermodynamical self-consistency which 
just expresses the demand that Legendre transformations between thermodynamical 
quantities, as they are derived from the partition function of the underlying theory, 
are not affected within the effective theory. In ,lj we have assumed a one-loop 
expression for the pressure to derive the evolution e(T). The purpose of this paper 
is to show that the one-loop evolution is exact for many practical purposes, that is, 
(thermal (quasi)particle) excitations in the electric phase are almost noninteracting 
throughout that phase^. 

The paper is organized as follows: In Sec.|21 we set up the real-time-formalism 
Feynman rules in unitary-Coulomb gauge and some notational conventions useful 
for organizing our calculations. In SecElwe sort out the diagrams that do contribute 
to the two-loop pressure for the SU(2) case. We discuss their general analytical form 
and defere hard-core analytical expressions to an appendix. Kinematical constraints 
on the off-shellness of quantum fluctuations as well as the center-of-mass energy en- 
tering a four-gauge bosons vertex are being set up and discussed. In the absence 
of external probes to the thermalized system these constraints derive from the ex- 
istence of a compositeness scale characterizing the thermodynamics of the ground 
state. In Sec.|3]we perform an analytical processing of the integrals associated with 

nonvanishing two-loop contributions to the pressure. In Sec. El we discuss the prob- 
^Some interesting physics does, however, take place shortly before the theory settles into its 

magnetic phase I . We discuss its implications for the large-angle 'anomaly' in the power spectrum 

of the cosmic microwave background in the last section of the present paper. 
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lems inherent to a numerical evaluation of loop integrals and their solutions. For 
the vacuum propagation integrals are either evaluated in a Euclidean rotated way 
and a subsequent imposition of the kinematical constraints or by performing e ^ 
limits numerically in the Minkowskian expressions. In Sec. E] we present our results 
graphically. In Sec.[7|we discuss and summarize our work and point towards its pos- 
sible phenomenological importance for the explanation of the large-angle anomaly 
observed in the CMB power spectrum jlOjj. 



2 Feynman rules and notational conventions 

In this section we set up prerequisites for our calculations. The two-loop diagrams 
for the thermodynamical pressure split into the contributions as displayed in Fig.^ 
There are local and non-local contributions. We will evaluate them within the real- 
time formalism of finite-temperature field theory JT]. For an SU(N) Yang-Mills 
theory the following rules apply: 

1. Each diagram is divided by a factor iV, where V denotes the number of ver- 
tices. 

2. Local diagrams are multiplied by a factor 1/8, nonlocal diagrams by 1/4. 

The three- and four-gauge-boson vertices ^[^^^bciPy ^' l) ^[4]abcd respectively 
(see Fig|21): 
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Figure 2: The vertices T^.^{,^^{p, k, q) and rf^'f^d- 

/ade/.ce(^7'^^/'^ - ^^'^^^7^'^)] (l) 

Since the effective theory has a stabihzed^, composite, and adjoint Higgs field 
characterizing its ground state, we shall work in unitary gauge where is diagonal 
and the pure-gauge background is zero (see [T| for a thorough discussion of the 
admissibility of this gauge condition). There is a residual gauge freedom for the 

unbroken abehan subgroup"^ U(l)^'\ A physical gauge choice is Coulomb gauge. 

■^The field (j) is shown to not fluctuate statistically and quantum mechanically [Tj. 
*This assumes maximal breaking by (f). 



4 



In unitary- Coulomb gauge each of the propagators for Tree-Level-Heavy/Massless 
(TLH/TLM) modes split into a vacuum and a thermal part as follows [T| ITT]: 



p2 _ ^2 



2n6{p'~m^)nB{\po/T\) 



(2) 



Dl^'Mip) 



D,u{p) 



-5abD^^{p) 



- + 2nS{p')ns{\po/T\) 



^IJ.U p2 



if fi = or u = 
else 



In Eq. (121) ub^x) = l/{e^ — 1) denotes the Bose- Einstein distribution function, and 
T is the temperature. We have neglected the propagation of the Aq field in the TLM 
propagator since we expect that this field is strongly screened - for T > 2 T^, where 
the TLH mass is sizably smaller than T, the Debye mass is ~ eT with e ~ 5.1, for 
T Tc there is an exponential suppression of this screening. 

With these rules at hand the two-loop correction to the pressure is given as 



AP — APi^^^i + AP^^^i^^^i 

o 4 



where the local contributions can be written as 



(3) 



1 fd^pd^k 



(27r 



|8 '^[4{abcd-^l^^,0'b{p)Dpa^cd{k) 



(4) 



For nonlocal diagrams we have 



1 [ d pd k X 

Yi J "(2^ ^[ZbM ^> -P - k)^Wsti-P^ -k,p+ k) X 



D\p^ar{p)D^,„^bs{k)Dy^^^t{-p - k) . 



(5) 



In Eqs.(0]) and D^^ab stands for both TLH- and TLM-propagators, and one has 
to sum over all combinations allowed by the vertices r[3^£jjj, and r|4]Qftcii" 

Let us now introduce a useful convention: Due to the split of propagators into 
vacuum and thermal contributions in Eq.Q combinations of thermal and vacuum 
contributions of TLH and TLM propagators arise in Eqs.(I3]) and We will con- 
sider these contributions separately and denote them by 

.pXYZ/XY .gv 

axf^Ylz/axl^Y ^ ' 

where capital roman letters take the values H or M, indicating the propagator type 
(TLH/TLM), and the associated small greek letters take the values v (vacuum) or 
t (thermal). 

3 Contributing diagrams for SU(2) 

In what follows we only investigate the case SU(2). It is clear that not all combina- 
tions of TLH- and TLM-propagators may contribute. This is due to the structure 
constants entering the vertices. For SU(2) they are fahc = ^ahc- As a consequence, 
the thirteen (naively) nonvanishing diagrams are 

^pHM ^ AP„^^ + AP„f^ + AP,f ^ + AP,f ^ 
APhhm = AP™^ + APj^r + AP^ + 

^pnUM ^ ^pUUM ^ ^pHUM _ (7) 



The number of allowed diagrams reduces further if one considers the strong coupling 
limit for the effective gauge coupling e (e > 0.5). This can be seen by virtue of the 
following compositeness constraint P]: 

|p2-m2|<|0|2 or pI + < \(f)\'^ (8) 

where the index E stands for the Euclidean rotated momentum. Eq. (jSI) expresses 
the fact that the ground-state physics is characterized by a scale set by |0| which 
determines the maximal hardness for the off-shellness of gauge-boson fluctuations 
the ground state can possibly generate. A Gaussian smearing of this constraint for 
Euclidean momenta introduces a ridiculously small effect since the variance for this 
distribution is, again, given by |</)p. Notice that only in unitary-Coulomb gauge, 
that is, the only physical gauge, it makes sense to impose the constraint (jH)). 

By the adjoint Higgs mechanism the (degenerate) mass of the two TLH modes 
is given as [T] 

m = 2e\(f)\ (9) 

where 

In Eq. (fTUj) A denotes the Yang-Mills scale. For later use we introduce a dimension- 
less temperature A as 

27rT , , 

A=— . (11) 

From the one-loop evolution of the effective gauge coupling it follows that e runs 
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into a logarithmic pole 

e(A)~-log(A-A,) (12) 

at Ac = 11.65 p. This is the point where the theory undergoes a 2""* order phase 
transition by the condensation of magnetic monopoles, and thus Ac corresponds to 
the lowest attainable temperature in the electric phase. 

We can scale out |0| in Eq. (jSJ. Then the Euclidean constraint becomes 

y^i^2qr(2i)2 < 1 (13) 

where = Since is always positive we conclude that only for e < 0.5 we 

do get a contribution from TLH vacuum fluctations in loop integrals. The plateau- 
value^ for e is, however, e ~ 5.1 as a result of the one-loop evolution P]. TLM 
vacuum modes do contribute, however, and we are left with the computation of 
AP^f ^, APjf ^, APjf ^ and AP^^^^ (AP^^^^ vanishes by momentum conserva- 
tion). 

There is one more kinematical constraint: For a thermalized system with no 

external probes applied to it, the center-of-mass energy flowing into a four- vertex 

must not be greater than the compositeness scale |0| of the effective theory. That 

is, the hot-spot generated within the vertex must not destroy the ground state of 

the system locally since the modes entering the vertex were generated by the very 
^This plateau indicates the conservation of isolated magnetic charge for monopoles contributing 

to the ground-state thermodynamics. It is an attractor of the (downward) evolution signalling the 

UV-IR decoupling property that follows from the renormalizability of the underlying theory. 
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same ground state elsewhere. This is expressed as 

|0p>|(p + fc)2| (14) 

where p and k are the momenta of the modes entering the vertex. As we shall see, 
Eq. ()14p leads to a strong restriction in the loop integration. 

To perform the contractions in Eqs. (jH) and (0) it is useful to exploit the transver- 
sality of the tensorial part D^y{q) of the TLM propagator from the start. The 
following four relations hold: 

D,u{q)q^ = 

DMa'"" = -2 



iq 



2 



D.M^k" = Pk-^^^' (15) 

The results for all relevant contractions are derived in the Appendix. 



4 Calculation of the integrals 

With the contractions of tensor structures at hand, we are now in a position to 
calculate all two-loop corrections. For AP^^^ this is done in detail, for the other 
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contributions we resort to a more compact presentation. We have 



2e^ / -^^-77- ( 24 - 6^ - 6— + 2 



{2tc)^ \ m? m? 

S{p' - m')nB{\po/T\) 6{e - m')nB{\ko/T\) . (16) 

In Eg. (I16|l both color indices a,c are summed over a, c = 1,2. The product of 
(5-functions can be rewritten as 

6{p^ - m^)6{k^ - m") ^ 



5(po + + m2)5(A;o - + m^) + 



The contraction V^^^^JI) ^y{j))D p„(}z) contains only even products of k and p (this 
is also true for the other contractions), like p^, k^ or pk. Thus, performing the 
zero-component integration over either 

5(po — a/p^ + m2) 5(A;o + ■\/k2 + m2) 
or 

(5(po + v^p^ + vn?) 6{ko - A/k2~+m2) , 

(signs in the argument of 5- functions opposite, crossterms) leads to the same result. 
This is also true for the two uncrossed products of 5-functions with equal signs. 
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After the integration is performed we may therefore set 



2 2 

p ^ m 



(pk) ±\/ p2 + Vk^ + m'^ — pk 



{pkf p^k^ + (p2 + k^)m^ + m^T 2pkv/p2 + mWk^ + m'^ + (pkf . 

(17) 

The upper case is obtained when the signs are equal, the lower case when they are 
opposite. 

Examining the integration constraint in Eq. (jl4p after the zero-component inte- 
gration over the products of 5-functions is performed shows that only the combina- 
tions with opposite signs must be evaluated: 

> \{p + ky\ = \pl-p^ + kl-k^ + 2poko-2pk\ 



> \2m^ ± 2v/p2 + mVk2 + m2 - 2pk| 



^ 1 > 2\{2ef ± ^/x^ + (2e) Vl/^ + (2e)^ - xycos^l 
where we have introduced (also for later use) dimensionless variables 

a; = |p|/|0|, i/=|k| 



z = cos^, x-'/' = ^. (19) 

27rT ^ ^ 

In Eq. (jl8|) 6 denotes the angle between p and k. We observe that for the "+" case 
the difference between the second and third term is always positive. And, because 
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of the first term, the whole expression is greater than unity in the strong couphng 
limit. Thus only the "— " case needs to be considered. This is also true for ^P^^' 
and AP^^J^'^ though the analytical expressions may look different. 



Applying the knowledge gathered in the Appendix, AP^^^ can be reduced to 

" J (27r)6 V rri^ J 

K( 2\x(i , /Tl^ 2\^B{\Pi)lT\)nB{\h/T\) 

S{po - + m^)6{ko + Vk2 + m^) . 

2v P + V k^ + 

Integrating over the zero components by using Eq. ()17|). we arrive at 



tt 



(2vr)6 2Jp^ + m^VWT 



12-2- 2 2 

After a change to polar coordinates and an evaluation of the angular integrals the 
remaining integration measure takes the form 2(27r)2|pp|kp(i|p|(i|k|(icos6'. As a 
last step we re-scale variables according to Eq. p9p . This re-casts the kinematic 
constraints of Eq. ()18|) into the following form: 

-l/2< (2e)2 - ^x^ + (2e) V?/^ + (2e)^ - xyz < +1/2 . (20) 

Our final result for AP^^ reads: 
TLH-TLH-thermal-thermal: 

^Pu = / dxdydz —====JL-==== X 



J + (2e) V?/' + (2f 

^W? W? 4(^V^ VI/ +(2 

riB (27r A-=^/2 + (2e)2)nB (27r A-^/^ ^^TW) 

12 



where the integration is subject to the contraint in Eq. ()20|) . The other two-loop 
corrections AP^^*^, AP^^ and AP^^^ are calculated in essentially the same way: 
TLH-TLM-thermal-thermal: 

We have 

nB{\polT\)nB{\k,/T\)5{p' -m^)5{e) 

where the sum is over a = 1, 2. Consider the integration constraint Eq. ()14|1 : 

W > \{p + qf\ = \pl-p^ + kl-k' + 2poko-2pk\ 
101^ > ± 2|k I a/p^ + m? — 2|p||k| cos^l 
^ 1 > \{2ef ± 2y^x^ + (2e)2 - 23;?/ cos ^| . (22) 

Again, the "+" case cannot be satisfied in the strong coupling limit (e > 0.5), so 
only the "— " case needs to be considered. Then AP^^'^ reduces to 
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AP™ = -2e' / -f-- -12 + 4^ + 2^ — x 

nB{\po/T\)nB{\k^/T\) 



2|k|^/p 



2 _i_ ^2 



5(po- Vp2 + m2)5(A;o+|k|) 



' (i|p| (i|k| d{cos6) -8 + 2 



(27r)4 y I I ^ ^: 

nB{Vp' + myT)nB{\k\/T) 



2e2T4 /• , , , of x^a - z^' 

dxdydzx y \ —8 + 2- 



(23) 



nB(27rA-3/Vx2 + {2ef)nB{2Tx\-^/^y) 
+ (2e)2 ' 

subject to the constraint Eq. (j221). 
TLH-TLM-thermal-vacuum: 
We have 

After the po-integration is performed the integration constraints Eqs. (jHl) and (fT^ 
read: 



\k'\ < 



b'-y'\<i 



(24) 



\{p + kf\< 



I (2e)^ + 7^ - ± 27 v/a;2 + (2e)2 - 2xi/ cos 0| < 1 



(25) 



Thus, the /cq- or 7-integration (7 is the re-scaled /co-component) cannot be performed 
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analytically. We have 



HM 
tv 



-2ie 



2 / d'^pd'^k 



riBWP^ + myT) 1 



, ^P^ sin 6^ 
12 + 4^7 + 2^^ — I X 



k^ 



^iPo - VP^ + 



-Aie 



pz, _|_ 

2 / (i|p| (i^o c^lk] (i(cos( 



p^k^ ( -8 + 2 



(27r)VP^ + 



p2 sin^ 6 



kl-k^ 



—Aie'^T^ f dx dy d'j dz 



xV ( -8 + 2^ 



(27r)5A6 J + {2ef 

{2TiX-''/^^x^ + {2ef)^ 



xHl-z^ 



(2e) 



riB 



^2 _ ^2 • 



(26) 



This is, however, not easy to evaluate numerically. To show the smallness of /S.P^ 



HM 

tv 



we resort to estimating an upper bound on the modulus of the integral in Eq. (j2f 
This is done by neglecting the center-of-mass energy constraint Eq. ()25p completely 
(but taking into account the constraint Eq. ((211)) and by integrating over the modulus 
of the integrand: 



tv 



d^pd^k 1 
(27r)7 P 

2 f d\p\ d{cos6) dk 



-o + Z I X 



e 



dx dz 



k 



P 



+ 2^ 



p2sin2^\ nB{\^/p^ + myT\) 



-8 + 2 \_ ,„ ^ I X 



(2e) 



X 



'p^^+rn? 



^x^ + (2e)2 



X 



nB(27rA-3/Va;' + (2e)2^ 



(27) 



In the second line of Eq. fl27|) /c has the meaning of /c = a/A; 
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TLH-TLH-TLM-thermal-thermal-vacuum: 

Here we have 



^pHHM _ _}_ f d pd kd q xt^,. , , .ypar , _ ^ 

^-^ttv - 2i J (27r)6 ^ [3]a63l^'''^'^M [3]a63l D 

Dxp{p)D^^{k)D,,{q) nBi\po/T\)nBi\ko/T\) x 
— m^)5{k'^ — m^)5{q + p + k) 



{k + py 

where the sum is over a,b = 1,2. Due to momentum conservation both kinematic 
constraints, Eqs. ^ and (fT^. are equivalent: 

|g|2 = + kf\ = + + 2pk\ < \^\^ 

|2(2e)^±2v/x2TT2ij2v/y2:r(2ij2-2a;ycose| < 1. 
This is the same as for AP^^^ , so only the "— " case needs to be considered: 

10p2 + lOA;^ + 16pk - 2^ - 2— - 



2 /> y-} 2 



AP 

"'^ 2z 7 (27r)6 



p'^{pk) k'^{pk) {pkY p^k^sin^^^ 

^ 2 ^ 2 2 ? 1 , ^ 

m"^ (p + kj^ 

10-34-34-84+4+4+ 

m"^ 

m'^ m'^ 
^^(a/p^ + m2/T)r;,B(Vk2 + m? /T) i 

X 



5{pq - v/p2 + m2)5(A;o + Vk^ + m2) . (28) 



Using Eq. (fTTjl . the part oc '^(p^,^)^^ in the square brackets after po and /cq integration 
reads 

p2 , k2 , y\.\l-z') , lp||k|;. 



12 + 4i47 + 4— + ^ ^ + 8^^^^x/p' + m2v/k^ 



m2 m2 
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where polar coordinates have already been introduced. 
The remaining part is 



-16 



p^k^(l 



+ 2 



IpI 



The propagator l/{p + k)"^ becomes after re-scaling 



(p + ky 2(2e)2 - 2y/x^ + (2e) V?/' + (2e)2 - 2xyz 



Thus, we have 



HHM 

ttv 



dx dy dz 



x'^y'^ 



X 



2A6 ; + (2e) Vl/' + (2e) 

nj,(27rA-3/2 7x2 + (2e)2) nB(27rA-=^/2^.^2 ^ (2e)2) 

(2e)2 - ^ (26)2^^2 ^ (2e)2 - xyz 
|l6p + ^2 ^ 2-^7x2 + (26)2 71/2 + (2e)2 + 



x''|/''(l -I- ^ 
(2e)2 

x2?/2(l 



-, 



+ 



(2e)2 

x2y2(^l _ 2;2) 

x2 + y2 _|_ 2x?/2; 



12 + 4 



X 



+ 4- 



(2e)2 (2e) 



+ 



+4 



(2e)^ 



xyz 



/x2 + (2e)2 7?/2 + (2e)2 



(29) 



5 Numerical integration 

The objective of this section is to numerically evaluate the expressions fl^ . (j^ . (P7j) . 
and (El. 

Two observations should already be pointed out here: 
(1) As it will turn out, ignoring the kinematical constraint Eq. ()25|) in the expression 
for AP^^^ gives an upper bound which is much smaller in modulus than the by- 



far dominating contribution subject to these constraints for A not too far above 
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Figure 3: Constraints on integration for AP^^: The lower bound in z is always — 1 
while the upper bound is not. White area: restriction is always satisfied, Zmax = 1, 
grey area: restriction can be satisfied, Zmax is given by z+ as in Eq. ()32|) . black area: 
restriction can never be satisfied. 

Ac. While for the former the exact implementation of the constraints is virtually 
impossible it is difficult but doable for the others. 

(2) The nonlocal correction has a singular integrand due to the TLM propagator 
being massless. 

Both problems are resolved in the following two sections. 
5.1 Constraints on integrations 

A straight-forward implementation of the kinematical constraints is to multiply the 
integrands with appropriate B-functions. This, however, cannot straight-forwardly 
be fed into a Mathematica program. Here, we demonstrate how the problem is 
tackled for AP^f^. 
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Eq. ()20|) can be rewritten as 

1 + 2(2e)2 - 2v/a;2 + (2e) Vl/' + (2e)' 



z = cos 9 < 



2xy 
and 



. = cos.>-^ + ^P->'-^V™^V^^^'^ (30) 

2xy 

where — 1 < -2 < 1. For the lower and upper integration hmit we therefore get 



lower limit: Zmm = niin[l, max [—1, y, e)]] 

upper limit: Zmax = max[-l, min [1, z+{x, y, e)]] (31) 

with the definitions 



-1 + 2(2e)2 - 2Jx^ + (2e) V?/' + (2e)2 

z_{x,y,e) = 

2xy 

Notice that Zmin{x,y, 5.1) always equals —1. A contour plot for z^{x,y, 5.1) is 
displayed in Fig. El This plot shows that the constraint hardly ever is satisfied. We 
observe that z+{x, y, 5.1) is smaller than —1 in the black area, greater than +1 in the 
white and inbetween these boundaries in the grey area. The integration is restricted 
to a small band around x = y only. Parameterizing this area leads to an upper and 
lower limit for the integration range iny = y{x) (depending on x). Looking at Fig.|2l 
one also sees that x runs from zero to infinity. For AP^J^^^ the constraints are the 
same. For AP^^*^ we have to re- adjust our definitions of the integration limits. The 
upper and lower limits of integration in z formally are defined as in Eq. (|31|) with 
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the difference that z± now are given as 

-l + 2(2e)2-2yVFTT2^ 



z_[x,y,e) 



2xy 



z+{x,y,e) = . (33) 

2xy 



5.2 Singular integrand in the nonlocal diagram 

For AP^J^^^ an additional problem arises. Consider the integrand: 



x^y^ nB{27cX-'^/^^x^ + (2e)2)nB(27rA-3/2^^2 + (2 



e 



^x2 + (2e)2^y2 + (2e)2 (2e)2 - ^x^ + (2e)2^y2 + (2e)2 - xyz 



(2e)2 (2e 



^2 + 7/2 + 2xyz 



Vx2 + (2e)2v/i/2 + (2e)2] } . (34) 



"(2e)2 (2e)2 (2e 

(2e) 

The first part in curly brackets has no singularity and can be integrated numer- 
ically without additional thinking. The part oc ^2^y2^2xyz '"^^ integrated 
numerically as it stands since it diverges at x = ?/ and z = —1. Complex analysis, 
that is, the residue theorem, can not be applied to this problem because we can not 
close the line integral at infinity due to the integration constraint. We therefore add 
ie (e > 0) to the inverse TLM propagator. One needs to prescribe a small value for 
e and check the numerical convergence of the integral in the limit e — 0. The results 
for A = 70, 200 are shown in Table l^^ The real part stabilizes while the imaginary 
part converges to zero. In our computations a value e = 10~^ is reasonable in view 
of available numerical precision. 
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e 


1 n — 6 
iU 


1 n-7 
iU 


1 n-8 
iO 


1 n-9 
iO 


1 X Ke i^rtty — i\j ) 


Z.iUzo 


z.iUoi 


Z.iUbo 


Z.iUoU 


T-^x Im AP™^(70) 


-3 X 10-^ 


-1 X IQ-^ 


-3 X 10-*^ 


-3 X 10-6 


T-^x Re AP,f/^^(200)(x - IQ-^) 


9.8744 


9.8850 


9.8882 


9.8892 


T-^x Im AP,f/^^(200) 


-1 X 10-5 


-5 X 10-6 


-1 X 10-"^ 


-4 X 10-^ 



Table 1: Numerical evaluation of the e-dependent part of AP^^^^^ (adding a term 
ie (e > 0) to the inverse TLM propagator in Minkowskian signature). Obviously, 
the real part of the integral is not sensitive to the value of e while the imaginary 
part tends to zero for e ^ 0. 

6 Results 

Having performed the numerical integrations, we now are in a position to present our 
results for each contributing diagram by plotting the ratio of two-loop to one-loop 
diagrams as a function of the dimensionless temperature Ac = 11.65 < lambda < 
200. Figs.m through [7| show the results. Notice that the one-loop result, see PP for a 
calculation, does not contain the contribution of the ground state. Notice also, that 
we kept e = 5.1 for all values of A thus ignoring the logarithmic blow-up of Eq. ()12j). 
Due to the exponential suppression for large e this yields an upper bound for the 
modulus of each diagram in the critical region. 
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Figure 6: An upper bound for the modulus of the ratio of ^AP^^'^ and Pj.ioop as a 
function of 11.65 < A < 200. 

Our computation indicates that the two-loop corrections are at most 0.2% of the 
one-loop result. The dominant contribution comes from the nonlocal diagram in 

Fig.in 



7 Summary and Outlook 

Our results can be summarized as follows: The picture of almost noninteracting 
thermal quasiparticles that was underlying the one-loop evolution of the effective 
coupling constant e in the electric phase of a thermalized SU(2) Yang-Mills theory is 
confirmed by the two-loop calculation of the thermodynamical pressure. The (tiny) 
modification of the one-loop evolution equation for e due to two-loop effects will 
be investigated in |12j. On a mesoscopic level this modification can be understood 
in terms of scattering processes off magnetic monopoles whose core size becomes 
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Figure 7: Ratio of iAP^f/^^ and Pi.i„„p as a function of 11.65 < A < 200. 
comparable to the typical wave length of a TLM mode for T \ Tc where 



r-p _ Ac A 



AcA _ 11.65A 
27r 27r 



denotes the critical temperature for the 2""* order transition to 



the magnetic phase. For T ^ Tc the magnetic charge of a monopole is too much 
smeared to be 'seen' by the TLM mode. This simple fact arises from the constancy 
of e for large temperatures and the core size or charge radius R{T) of a monopole 
being approximately its inverse mass M pQ 

h-rrT 

(35) 



R{T) ~ M-\T) ~ e 



A3 



Thus the quick die-off of the two-loop correction to the pressure at large T (compare 
with Figs. m through [7| and ignore the fact that our estimate for AP^^*^, see Fig.l^l 
is too rough for large T due to the omission of the vertex constraint and that a 
small infrared effect survives for large T in Fig. El due to the masslessness of the 
TLM mode). The mechanical analogon for this situation is as follows: Imagine 
a box filled with heavy lead balls being at rest and light ping-pong balls moving 
around them. Now, switch on an interaction between the two species (wavelength 
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of TLM mode becomes comparable to charge radius of monopole for T \ T^). 
This will thermalize the system. However, the average momentum that is deprived 
from the ping-pong balls and added to the lead balls does not have an effect on 
the partial thermodynamical pressure of the latter since their momenta only probe 
the exponential tail of their Bose distribution. On the other hand, a decrease of the 
average ping-pong-ball momentum sizeably decreases their partial thermodynamical 
pressure. This is seen in Fig.EJby the (negative!) dip of the dominating two-loop 
correction. 

Despite the large value of e ~ 5.1 the smallness of two-loop corrections emerges 
from the existence of compositeness constraints which in turn are derived from the 
existence of a nontrivial ground state. We expect no major complications when 
generalizing our computation to SU(N). The situation is somewhat reminiscent of 
M = 2 supersymmetric Yang-Mills theory where the perturbative (3 function for 
the gauge coupling is exact at one loop JH]- The important conceptual difference 
is that the one-loop exactness in the supersymmetric case is inforced by a strong 
symmetry while in our approach to the A/" = Yang-Mills theory the identification 
of the essential degrees of freedom makes the interactions thereof almost vanish. 
We expect that the loop expansion of the thermodynamical pressure of an SU(N) 
Yang-Mills theory is not asymptotic but converges very quickly. 

An important application of our results arises: If the photon is generated by 
an SU(2) Yang-Mills theory of Yang-Mills scale A ~ Tcmb ~ 10~^eV being at the 
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boundary between the magnetic and electric phases but on the magnetic side^ then 
hght, being released at the time of decoupling of the CMB (deep within the electric 
phase of SU(2) 

cmb)) must have travelled through a 'lattice' of scattering centers 
{dual magnetic, that is, electrically charged monopoles) shortly before the Universe 
settled into the CMB dip where the monopoles are condensed into a classical field 
pp. This effect is seen in Fig. I?] by a decrease of the dominating two-loop correction 
to the pressure for T approaching Tc (that is Tcmb) from above. The observable effect 
should be a cosmic Laue diagram with a large quadrupole contribution and manifest 
itself in terms of a large-angle 'anomaly' in the power spectrum of temperature 
fluctuations in the cosmic microwave background. Such an 'anomaly' indeed has 
been reported by the WMAP collaboration [TUj . 
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^Only there is the photon precisely massless and completely unscreened: a situation which is 

dynamically stabilized by a dip of the energy density at Tcmb CO- 
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Appendix 

Here we evaluate the contractions of the tensor structures as they appear in Eqs. 
and ©• Exploiting Eq. (jl5p . the contractions for local contributions are: 
(1) Local, TLH-TLH: 



PfiPu \ f kpk, 



9fJ.l^ 9 j ( 9p<J 9 



-2^e^(24-64-64 + 2^-2^) • (36) 



(2) Local, TLH-TLM: 



eaceeMe{9'''9"' ' 9^"" 9""') + ta,,tU9''''9''' - 9^'9n] x 
'eaceeace [29^" 9"" ' 9'" 9"' - 9'' 9'"'] X 



{9pu-^)Dp^{k) 



_2.eM-12 + 4^, + 2P!^~) . (37) 



27 



In Eq. fj37|) 9 denotes the angle between p and k. For the nonlocal diagram we 
obtain: 

= ehabcerst [g^'^ip - ky + g^\k - qf + g^\q - pY] x 
[gP<r^p - ky + g'^^{k - qy + g-^P^q - py] x 

5ar5bs5ct {gxp - ^) {g^.. - ^) D,,{q) . (38) 

For not loosing track, we split the calculation into terms oc e^, oc ^ and oc ^ and 
keep D uncontracted in a first step. The contraction of structure constants eabz^ah'i 
gives an additional factor 2. 
Term oc 2e^: 

^-(A; _ qy + ^-P(g - py]gx,g,aD,r{q) 
= [gpaip - ky + gl{k - q)p + (7;(g - p),] x 

_ ky + <7'^^(fc - qy + <7-''(g - p)'^] D^M 

= [A{p - ky{p - ky + {p- ky{k - qy + {p- ky{q - py + 

(k - qy{p - ky + {k- qfg^^ + (g - py{k - qy + 
(g - Pfip - ky + {k- qy{q - py + (g - pfg^^] D^M 
= [2pV + '^k'k'' - Qp-'k^ + (g - pfg""^ + (k - qfg""^] D,^{q) 

2{q-pf-2{k-qf. (39) 
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Terms proportional to q" or have been omitted after the second-last equal sign in 
Eq. because, when contracted with D^T-{q), they vanish. Again, using Eq. (fT^ 
the expression after the last equal sign in Eq. (jH^ easily follows. 
Next we look at the two terms proportional to 2^ (compare with Eq. (j38j) ). 
The first one is: 

[g^^ip - ky + g^'^ik - q)' + g^\q - p)^] [g^-{p - kf + 
^-(fc - qf + g-f{q - pY]g^pk,k,D,M 
= [kpk,{p - ky + rk,{k - q)p + k{q - p)g';k,] X 
[g^"{p - ky + g"^{k - qy + g^^{q - py]D,r{q) 

= [k\p - ky{p - ky + k{k -q){p~ kyk^ + k{q -p){p- kyk^- + 

k{k - q)k''{p - ky + {k- qfk^k-" + k{q - p)k''{k - qy + 
k{q - p)k^{p - ky + k{q -pyk- qyk^ + [k{q - p)fg'"'] D^M 
= [Pp'p'- + q^k'k"- - 2{kp)p'k^ + [k{q - p)fg''^] D,r{q) 

-2[kiq-p)]'. (40) 
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The second term oc 2^ either is obtained by a direct calculation or by just exchang- 
ing p ^ k in Eq. (PUI) : 

[g^^ip - ky + g^'{k - qf + g''\q - pY] [g^^^p - + 
^-(A; - q)P + _ pY]g,„p^p,D,M 

2\p{q-k)f. (41) 

2 

Finally, the term oc 2^ is given by 

_ ky + g^''{k - qf + g''\q - pf] [g^-^ip - kf + 
^-(A; - q)P + ^^^(g -p)-]pAPpA:;.fc.5..(g) 
= [(pfc)Pp^f7(p - ky + p(A; - q)ppk„k'' + k{q - p)ppk„p''] x 
[g^^ip - A;)- + g'^^ik - qY + ^?-''(g - py]DUq) 

= [{pkyp - ky + \p{k - g)]r + \k{q-p)Y\ 

{i^kyp - A;)" + \p{k - g)]F + \k{q-p)Y\DuM 
= [{kq)Yp^ + {pq^k^k^ - 2{pq){kq)p''k^]D,^{q) 

- (*^)^(p^-^) + (M)'(k^-^)- 

2(M)M(pk-<^«)^ (42) 
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Now, adding up Eqs. ljH^ through (jl^ (taking care of the correct signs), we have 



Eq. (jSHD = [Eq. (jSHl) - Eq. (gOD - Eq. (jUD + Eq. 



{kqf 



nv 



P 



pY-2{k-qy + 



[2 



(pqr 

(kq)2 



(43) 



Using momentum conservation at the vertices, that is g = —p — k, we find: 



P 



(kq)^ 
|qP 



pk 



(pq)(kq)\ p^k^sin^6' 
|qp / (p + k)2 



(44) 



And thus, 

(3) Nonlocal, TLH-TLH-TLM: 



2e' 



10/ + loe + i6pk - 2^ - 2^ - "'''^^^^ 



16 



(pk) 



2 , in;.2 _^ 

2 ^2l,2„;„2 



4 



EJii|L!(io-3i^-3i^-84 + i^ + 

(p + /cj^ V m"^ 



P^{pk) .k'^ipk) ApkY rtP^k"^ 



+ 4 



+ 4 



+ 4- 



+ 2^ 



(45) 



Here, we have used the fact that T{—p, —k, —q) = — r(p, k, q). 
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